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On the Approximate Arithmetical Solution by Finite Differences 
of Physical Problems involving Differential Equations^ ivith 
an App)lication to the Stresses in a Masonry Dam'. 

By L. r. ElCHARDSON. 

(Communicated by Dr. E. T. Glazebrook, F.E.S. Eeceived in revised form 

November 2, 1909, — Eead January 13, 1910.) 

(Abstract.) 

In order to deal with irregular boundaries, analysis is replaced by 
arithmetic, continu.ou.s fu.nctions are represented by tables of numbers, 
differentials by central differences. Problems then fall into two classes : — 

(A) The relation between the equation obtaining throughout the body, and 
the boundary condition is such that the integral can be stepped out from a 
boundary. This class includes equations of all orders and degrees. It has 
been treated by arithmetical differences by Eunge, W. F. Sheppard, Karl 
Heun, W. Kutta, and Eichard Ganz. Examples of a specially simple method 
are given. 

(B) The integral must be determined with reference to the boundary as a 
whole, as in Dirichlet's problem. The method given has only been worked 
out for a limited group of linear equations, — namely for those in connection 
with which a function analogous to potential energy exists, which is a complete 
minimum when, and only when, the difference equations are satisfied. Under 
this condition the difference between the finite-difference-integral (j)^ and a 
function cj^i of the independents, having the correct boundary conditions but 
otherwise arbitrary, can be expanded in the form cj^i — ^u = ^AkPk, where 
the Ai...A^ are constants and Pi...P^ are "principal modes of oscillation" 
defined by WPk = ^/c^P/o where Wcj^u = is the difference equation to be 
integrated and X^ is a constant. Now we start with the table of numbers <^i 
and calculate 39'<^i. Then as lB'(f)u = 0, we have 

Multiplying both sides by some number ai~^ and subtracting from ^i, and 
altering the boundary numbers so that the boimdary condition is still satisfied, 
we have a new table which may be called ^2; and 

^2 — (?^w = ^Aa- (1 — oii~^Xk^) Tk. 

Eepeating the process with a2...a?n, we get 

^m+l — ^w =XAk (1 — ar^A//) (l--a2~^XA:^)...(l-~am""%^)PA. 
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Now a function I exists such that ^IP/ = 1, ^IPzP/^ = 0, where S? 
denotes a summation throughout the region. Therefore 

ISTow by a sufficient number of suitably chosen ce's the polynomial in \^ on 
the right can be made small throughout the range from Xi^ to X^^. Therefore 
the error of (pm+i can be made small ; for, since I is one-signed, it is measured 
by the L.H.S. The process is arithmetical. 

Under certain conditions the error due to finite central differences is of the 
form e2li^ -\- eji^ -\- bqIi^ -\- Qto., where h is the co-ordinate difference and the a's 
are functions of position independent of h. If the integral has been found 
for two or more sizes of h, more exact values of it can be extrapolated by 
this formula. 

These methods have been applied in the paper to calculate the stress- 
function in. a masonry dam. 



The Initial Accelerated Motion of Electrijied Systems of Finite 
Extent^ and the Reaction prod%tced hy the Resulting Radiation, 

By George W. Walker, A.E.C.Sc, M.A. 

(Communicated by Prof. Sir J. Larmor, Sec. E.S. Eeceived December 17, 1909, 

— -Read January 27, 1910.) 

(Abstract.) 

When the exact equations of motion of a dynamical system are known, it 
is in general possible, by a well-known process, to determine the initial mode 
of change from a steady to a variable state, even when the primary equations 
cannot be completely integrated in the general case. The primary object of 
this investigation was to show that the same process is applicable to the 
equations of motion of finite electrified systems. That the results would be 
applicable to the question of the electric inertia of " electrons " was constantly 
kept in view ; and it was felt to be undesirable to make any approximation 
depending on extreme smallness of the electrified system. 

The various expressions hitherto used for the electric inertia of " electrons," 
at speeds comparahle with that of radiation, depend, I think, without exception, 
on consideration of the energy of a steady state. The " quasi-stationary " 
principle assumes that, when the energy in a steady state is known, it is 



